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Abstract 

In case of non-dissipative tissue the inverse problem of thermoa- 
coustic imaging basically consists of two inverse problems. First, a 
function depending on the electromagnetic absorption function, is 
estimated from one of three types of projections (spherical, circular 
or planar) and secondly, the electromagnetic absorption function is 
estimated from (f>. In case of dissipative tissue, it is no longer possible 
to calculate explicitly the projection of from the respective pressure 
data (measured by point, planar or line detectors). The goal of this 
paper is to derive for each of the three types of pressure data, an 
integral equation that allows estimating the respective projection of 
(p. The advantage of this approach is that all known reconstruction 
formulas for <f> from the respective projection can be exploited. 

1 Introduction 

The goal of thermoacoustic imaging is to estimate the electromagnetic ab- 
sorption function of soft tissue so that the strong contrast between the 
electromagnetic absorption of cancer and normal tissue can be exploited 
(cf. [T3l [151 [HI ED EH EI])- Another advantage of thermoacoustic imag- 
ing is that the technological advance allows measuring different types of data 
(cf. [251 El tD EI]) and therefore various explicit reconstruction formulas can 
be applied (cf . [51 El EH E51 E3 EES ED E3] ) • In order to increase the resolution 
of thermoacoustic imaging, we derive integral equation models that take dis- 
sipation into account (cf. [2J,[T9lED]) and allow the application of all known 



reconstruction formulas. In the following we explain this more precisely in 
case of point detectors. 

First we discuss shortly the direct problem. The direct problem mod- 
els the propagation of a pressure wave in dissipative tissue generated by 
very short heating of the tissue due to a laser. In case of homogeneous 
and isotropic dissipation, the direct problem can be modeled as follows (cf. 
Theorem 4.2.1 in [8]) 

(1) p = G * x>i / (* Xjt space-time-convolution) , 

where G corresponds to an attenuated spherical wave with origin in (x, t) = 
(0, 0) and / corresponds to the source term of p. In the absence of atten- 
uation, G is the retarted Green function of the standard wave with source 
term /. If G is modeled appropriately, then / is the same sorce term as in 
the absence of dissipation (cf. Section [2] and [2~T]). i.e. 

(2) /(x, t) = *—~q, - <K X ) (I short time impulse) , 

where depends on the electromagnetic absorption function, say = F(a e i). 
The essential requirement for G is causality, i.e. the speed of the wave front 

(cf. H21) 

{(x, T(x)) G R 4 | T(x) := sup{t > | G(x, r) = for all r < t} 

has to be bounded from above. Here 7X|x|) denotes the travel time of the 
wave front. Below we will see that causality plays an important role in 
thermoacoustic imaging. 

The respective inverse problem can be formulated as follows. Let V denote 
a surface covered with point detectors surrounding the tissue. From (pP) 
and (j2]), we obtain the integral equation 

(3) /G(x-y,t)0(y)dy = p(x,t) for x G T, iG[0,T], 

where G := G *t and T is sufficiently large. The inverse problem cor- 
responds to the estimation of in integral equation ([3]) and the subsequent 
inversion of = F(a e i). 

In this paper we show that integral equation can be reformulated to 

oo 

(4) J iV(^,^ / ) J Rs P (0)(x,t / )dt , = />(x,t) (iVGC°°(R 2 \{(o,o)}), 
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where R sp {<p) denotes the spherical projection of 0, i.e. 

(5) i? sp (0)(x,t) := J 0(x')dA 2 (x') x£l 3 ,t>0. 

8B t (x) 

Here dA 2 denotes the Lebesgue measure on M. 2 . The causality of G ensures 
that N(t, t') — if t' > t, i.e. the upper integration limit in (j3J) can be 
replaced by t. This means that the pressure function does not depend on 
information from the future. For some surfaces T, can be reconstructed 
from R sp ((j)) via exact or approximate explicit reconstruction formulas. For 
example, if T is a sphere, then (cf. e.g. [T^ [To]) 



div / nfx' 



(6) W = 

s 2 



0_ 

ot 



(^ p (0)(x',t)) 



dA 2 (x' 

t=|x'-x| 



where n denotes the exterior normal of S 2 . In the absence of attenuation, 
if I(t) = 5(t), then the projection R sp (4>) can be calculated explicitly from 

(cf. nm) 

._. d /i? s „(0)(x, t)\ 

(7) — I — — =po( x , t) (po unattenuated data) . 



dt V 4vrt 

However, in the presence of attenuation, a set of Id— integral equations 
like (jl]) has to be solved to obtain the projection R sp ((j)). Moreover, we see 
that if an efficient and accurate explicit reconstruction formula for from 
data R sp ((f>) exists, then the iterative estimation of via integral equation (J3j) 
is unfavorable. The goal of this paper is to derive Id— integral equations that 
relate spherical, circular and planar projections to attenuated pressure data 
measured by point, line and planar detectors, respectively. We note that this 
is a theoretical paper that is not concerned with numerical experiments. 



The outline of this paper is as follows. The direct problem of thermoa- 
coustic imaging of dissipative tissue is modeled in Section |2j In Section [3] 
an integral equation model is derived that allows the estimation of the unat- 
tenuated pressure data in case I(t) = S(t). Then, from this model, the 
basic theorem of thermoacoustic imaging of dissipative tissue is derived in 
case I(t) 7^ 5(t). Subsequently, in Section HJ integral equation models for 
the three types of projections are derived that allow applying the explicit 
reconstruction formulas of thermoacoustic tomography. In the appendix it 
is proven that the wave equation modeled in Section [2] has a unique Green 
function satisfying the "initial condition" G|t<o = and causality, defined as 
in Section EJ 
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2 Modeling of the direct problem 



In this section we model the direct problem of thermoacoustic imaging of 
dissipative tissue. First we summarize essential facts about wave attenu- 
ation and causality and then model the stress tensor and the source term 
for thermoacoustic imaging of dissipative tissue. For more details we refer 

to [121 HZl HS1 E31 12D- 



Causal wave attenuation in tissue 

Consider a viscous medium that is homogeneous with respect to density, 
compressibility and attenuation. Then pressure waves obeying a complex 
attenuation law a* = a* (a;) satisfy the wave equation (cf. [T2] ) 

(8) V 2 p 7 (x, t) " + ^ f ) Pj(x, t) = -/(x, t) , 

where Co > is a constant and is a time convolution operator with kernel 
defined by 




(" denotes Fourier transform) . 



The functions a* : R — > C is called complex attenuation law and a := Re(a*) 
is called (real) attenuation law. In order that p 1 is real valued a must be even 
and Im(a*) must be odd. Attenuation occurs only if a is positive. Moreover, 
the function a* is restrained by the requirement of a bounded wave front 
speed of the Green function G 7 of (JS} (causality). If the wave front speed 
this is bounded from above by C\ 6 (0, oo), then this is equivalent to 

G 7 (x,t + T(|x|)) = if t<0 forallx^O, 

where T(|x|) G (|x|/ci, oo) is the travel time of the wave front from point 
to point x. In case of a constant wave front speed Cq, causality is equivalent 
to 

(9) K(x,t) := 4tt|x|G 7 ^x,t + ^ = if t<0 for all x ^ . 

since T(|x|) = |x|/c . We note that Kfcu) = e~ a * (w) ^/V2tt (cf. Theo- 
rem [9] in the appendix). 
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Remark 1. In the literature often a less strong definition of causality is used 
that demands the existence of a (retarted) Green function that vanishes for 
t < 0. However, this requirement is not related to the speed of the wave front. 



According to experiments the real attenuation law of a variety of viscous 
media similar to tissue satisfy (at least approximately) a frequency power 
law (cf. PllEZj). 

a(cj) = a |w| 7 for 76 (1,2], a = const. 

This led to the complex attenuation laws (cf. [221 [231 [2H [26] ) 



(10) a„(u) = a (-iw) 7 with a 



cos(7T7/2) 



which violate causality for the range 7 G (1,2] relevant for thermoacoustic 
imaging (cf. appendix). For thermoacoustic imaging we propose the follow- 
ing complex attenuation laws 

/ s dn (— i id) , 

(11) «.(") = 7= . = j ( T e(l,2,r >0 , 

C VI + (-1T W) 7 " 1 

where the square root of 1 + (— irow) 7-1 is understood as the root that 
guarantees a nonnegative real part of a*. 

Remark 2. We note that the above complex power function is defined by 

(12) w 7 = e "/( l °9( r )+^) f or w = r j<P e c~ , 
where C~ := C\{z E C | Re{z) < 0, Im{z) = 0}. 

Remark 3. In the appendix we prove that the Green function G 7 of |5]) 
with a* defined as in ( T7D|) does not satisfy the causality requirement |PP, 
but if a* is defined as in / T77]) then (EP is satisfied (cf. [12]). Models / T77]) 
are appropriate for thermoacoustic imaging, since causality is satisfied and 
Re(a*) is approximately a frequency power law with exponent 7 G (1,2] for 
small frequencies (cf. Figure\V$. If 7 = 2 then Re(a*) is the thermo-viscous 
attenuation law. 



Remark 4. All mathematical results in this paper are valid if a* defined as 
in (Q7p is replaced by a complex attenuation law that satisfies (Tj|) and has a 
monotonic increasing real part. 
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Figure 1: Comparison of the real part of (fTTi) with «o := \ C os(n^%)\ (dashed 
line) and the power law a(u) = |t u;| 7 for 7 = 1.5. For liquids: r = 
10~ 6 MHz (left picture) and for gases: To = 10~ 4 MHz (right picture) 
(cf. [llj). Experimental demonstrations of the power law are performed for 
the range - 60 MHz (cf. e.g. [23]). 

Modeling of the stress tensor and the source term 

Since equation (JSJ is an integro-differential equation, it is not self-evident 
that the source term / has the same structure as in the absence of wave at- 
tenuation. In order to model the source term we first model the temperature 
dependent stress tensor. 

According to [16] the stress tensor of an invisid fluid can be modeled by 

= ~ c o Po a th (T - T ) Sij - pi Sij where p 1 = p + c\ (p - p ) • 

Here po, po, cq, T are reference values of the pressure, the density, the sound 
speed and the temperature, respectively, and a t h is the thermal expansion 
coefficient. Wave equation (jHJ) shows that a dissipative medium has a long 
memory (cf. [3]) and therefore we claim that attenuation is accounted for if 
o~ij is replaced by Ki * t a it j, where * t denotes time-convolution and K 1 is an 
appropriate kernel. We show that this ansatz implies wave equation (jHJ) with 
the same thermoacoustic source term as in the absence of wave attenuation. 
(This calculation can be considered as another justification of model equa- 
tion Since = cr 2 ,2 = cr 3)3 correspond to the negative total pressure 
in the absence of dissipation, our claim implies 

(13) p = po + Cq Ki * t (p - po) + Cg po a t h K\ *t (T - T ) 

with po := po J R K x (t) dt = V2^Po Ki(0+). 
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If the activation time of the laser is very short, say At, then the absorbed 
heat content of Q at t = At can be modeled by (cf. [B]) 



AQ(x, At) = a el e (x) / e"^ I™' dx s I(t) At 
=: a d I (a el ,x) I(t) At, 

where ct e \, eo and x s denote the electromagnetic absorption coefficient, the ab- 
sorbed energy flow density and the origin of an electromagnetic point source, 
respectively. Moreover, / = I(t) denotes a short time pulse. On the other 
hand the heat content of the mass within Q due to the rise of temperature 
is 

AQ(x, At) = m n c p (T(x, At) - T ), 

where c p denotes the specific heat capacity of at constant pressure. In 
summary we get 

dt m n c p 

Since a state equation refers always to the local comoving system, the total 
time derivative d/dt appears. From (|T3i) we infer the state equation 

. . 1 dp T ^ dp T ^ dT 

(15) ^ = Kl * t dt +Po(XthKl * t dt- 

If the sound speed is not to large, then d/dt can be replaced by d/dt. In the 
following we assume that this simplification is appropriate. The linearized 
equation of motion and the linearized continuity equation 

<9v dp _ 

Po ^7 = - Vp and — + p V ■ v = 
dt dt 



imply 

9t 2 



From this, f JT5|) and f|14jl we infer 



(16) i^(t) V 2 p(x,t) - ~ 2 = -JTi(f) * t /(x,t) 



with 

/= a,^K-) a[ volume of ^) . 



c p dt 
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Comparison of (jSJ) and ffTSj) in Fourier space, shows that both equations are 
equivalent if 

(17) ^^-Lf^-V, fc(u):=io,(u) + -. 

V2 7T \c k(uj)J c 

Now we can formulate the direct problem. 



The direct problem 

Let T,Ti > with T x « T and e.g. G L 2 (fi) and J G £(R) with 
supp(J) C [0, Ti]. The direct problem of thermoacoustic tomography is to 
solve wave equation f)16p on M 3 x (0, T) with 

(18) /(x,t):=0(x)/'(t) where := I 



such that 



(19) P 7 l t <o = °- 

Let S(t) denote the delta distribution. Frequently, it is assumed that I(t) = 
5{t), since the speed of light is much larger than the speed of sound. Then 
we have 

(20) /(x,t) = 0(x)5'(t) 
and corresponds to an initial value function. 

Remark 5. In the appendix we show that wave equation [To]) has a unique 
Green function G 7 f/(x, t) = <5(x) 8(t) ) satisfying G 7 |t<o = 0. 



3 Derivation of the basic theorem 

Before we explain the purpose of this section we introduce some important 
notions and assumptions. 

(Al) Without loss of generality we assume c = 1. a* is defined as in (ITTj) if 
7 G (1, 2] and a* := if 7 = 0. Moreover, k is defined as in ( TTT|) . 

(A2) p 7 is the unique solution of (IT5|) with f lT9"j) and /(x, t) = </>(x) I'(t) and 
po is the unique solution of ( IT6|) with ( fl9|) and /(x, t) = 0(x) 5'(^)- 

If i"(t) 7^ 5(t), then we assume that the signal satisfies 



8 



(A3) I e C(R) and supp(J) C [0,Ti] for some T : > 0. 



If I(t) = 5(t) then p Q = p . The uniqueness of the Green function of ffTB]) is 
proven in Theorem [H] in the appendix. 

The goal of this section is to derive an integral equation that allows the 
estimation of p from the data p T This will be the basic theorem of ther- 
moacoustic imaging of dissipative tissue. 

The following two lemmas provide important properties of a distribution 
that plays an important role in deriving the basic theorem. 



Lemma 1. Let (Al) be satisfied and 

1 u 



(21) 



M 7 (u,t') :-- 



/27T k(u) 



ik(uj) \t'\ 



(u, t' e 



a) M 7 satisfies for t' > 
(22) 



,. i)-M. i) 2 \L 
Ki * t —^P- ^-x 2 = with 



dt' 2 



dt 2 



dt' 



-5'{t) 



t'=0 



b) lfa*^0 then M 7 G C°°(M 2 \{(0, 0)}). 

c) I/a* = then M 7 (t,t') = 5(t - \t'\). 

Proof, a) In Fourier space, equation (122]) and definition (12TI) read as follows: 
(-iio) 2 d 2 MJuA') 



'7n 



dt' 2 



-iw) 2 M 7 (u;,^ , ) = 



= (i k(u)) 2 M 7 (w, If) t'>0. 

From both identities, it follows that M 7 (i, if) satisfies the first identity in (T221) . 
The second identity in (1221) follows from (I2ip : 



dM y (t,t') 



dt' 



T- 1 



t'=0 



1W 



x 2 7r 



(t) = -5 f (t) . 



b) Let a* 7^ 0. For fixed \t'\ = T > 0, M 7 can be considered as oscillatory 
integral (cf. Section 7.8 in [8]) 



a(t,u) u(t)dtdu u e Cg°(R) 



with phase function (j)(t,u) = — tu in T = Ix (M\{0}) and a(t,u) = 
"7§^ TJu) e lfc( ^ T , since a is C°° and rapidly decreasing with respect to oj. 
From Theorem 7.8.3 in [8 J and 

^-(t,u) = iff t = 0, 

it follows that singsupp (I^a) ^ {0}. That is to say M 7 (£, T) is C°° for t > 0. 
Similarly, if i = T is fixed'then it follows that M 7 (T,t') is C°° for > 0. 
This shows that M 7 is C°° on M 2 \{(0,0)}. 

c) The claim follows from k(u) = u (a* = 0) and property ^{dlt — t')}(u) = 
e lu)t j\p2/K. This concludes the proof. □ 

Lemma 2. Let (Al) be satisfied with a* ^ and M 7 6e defined as in l[21\) . 
a) lft' = then 

M 7 (t, 0) = for t < . 

6j Ift' 7^ £aen t i— > M 7 (t, £') a rapidly decreasing C°°— function on K. and 
saizs/ies 

(23) M 7 (t,t') = /or t<\t'\. 

Proof, a) From (ED with t' = and (TTT]) we get 



V^M>, 0) = " = V /1 + ( - iT0U;)7 " 1 (ao > 0) . 

Here the power functions are defined on C\(— oo, 0]. Since a +\/l + ( — i r o 2 ) 7 ~ 1 
maps 

{zGC| Im(z) > 0} to {wGC| Re(w) > a } , 

we see that M 7 (z, 0) is holomorphic for all Im(z) > 0. Hence conditions (CI) 
and (C2) of Theorem [6] in the appendix are satisfied. Moreover, we see that 
there exists a polynomial P such that 

\M 7 (z,0)\ < P(\z\) for lm(z)>0. 

According to Theorem in the appendix, M 7 (t, 0) vanishes for t < 0. 
b) Let a* ^ 0. If t' ^ is fixed, then M 7 (u;,t') is rapidly decreasing and 
C°°— and thus M 7 (i, t') is a rapidly decreasing C°°— function on R. Hence 
property ( 123]) is satsified if 

(24) M 7 (t+ |t'|,t') = for t<0. 
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Because of (l2~Tj) and the second definition in (|T7|) we have 

M 7 (t + \t'\,t') = 7- 1 {M 7 ( W ,0 ^ '} (*) = ^ {^y } (*) 

According to part a) of this proof and Theorem [8] in the appendix both 
functions 

vanish for t < and therefore their convolution vanishes for t < 0, too. This 
proves property f[24|) and concludes the proof. □ 



Let n, T2 G R. We recall that if fi and /2 are two distributions with support 
in [n, oo) and [r 2 , oo), respectively, then /i * f 2 is well-defined and (cf. [8]) 

(25) supp(/i * / 2 ) C supp(/i) + supp(/ 2 ) C [n + r 2 , oo) . 

Now we derive the first integral relation that relates the attenuated pressure 
data with the unattenuated pressure data in case I(t) = S(t). Note that in 
this case p = p (cf. (A2)). 

Theorem 1. Let (Al), (A2) with I(t) = 5(t) be satisfied and M 7 be defined 
as in (dip. Then p solves the integral equation 

oo 

(26) p 7 (x, t) = J M 7 {t, t') po(x, dif . 

o 

i/ere t/ie upper limit of integration can be replaced by t. 

Proof. Let t! > 0. From properties (BIZ) . ([25]) and (|Z5)» . it follows that 

(27) Jj^t) * t - ^ tJ =0 if t < H and n G {0, 1} . 
Moreover, p satisfies (7 = 0, K\{t) = 6(t)) 

V%(x,0 = ^^ fort'>0, 

(28) ^ 

Po(-,O+) = and J£(.,Q+) = 0. 



9* 



For convenience we introduce the notions 



(29) Dp := Ki(t) * t (V 2 p) - |f 
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and 



oo 

g(x,t) := J M^t,t')p (x,t')dt' . 



Kernel property (1231) implies that the upper limit can be replaced by t. 
From (|22|) we get 



oo 

□g(x,t) = J K x {t) * t 



M 7 (t,f) V%(x,t') - Po(x,t') 



dt' 



which simplifies with integration by parts together with fl2T]) and fl28|) to 

0M 7 (t,t') 



(30) ng(x,t) 



dt' 



(x) 5(t) * t ^(t). 



J t'=0 



Because p 7 is the unique solution of (fi6i) satisfying (fl9l) (cf. Theorem [9] in 
the appendix), it follows that p 7 = g. This concludes the proof. □ 



From Theorem [T], assumption (A3) and f[2o^) . we get: 

Theorem 2. Let (A1)-(A3) be satisfied and M 7 be defined as in l[21\) . Then 
po solves the integral equation 



(31) 



oo 

p 7 (x, t) = y [/(t) * t M 7 (t, O] Po(x, dt' . 



The upper limit of integration can be replaced by t. 



4 Integral equation models for projections of 

In this section we derive integral equation models for the three types of pro- 
jections in thermoacoustic tomography. For various situations, the function 
<p can be calculated via explicit reconstruction formulas from projections of </> 
(cf . El [H H5J EHl EH EH [33] ) . In case of reconstructions with limited- view, 
we refer to [32J and the literature cited there. 

In this section we use the following notions and assumptions: 

1) The region of interest (tissue) 11 C M 3 is a subset of the open ball 
.Br (0) with radius Rq > 0. 

2) dA m denotes the Lebesgue measure on M. m for m £ {1, 2, 3}. 
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Case of Point detectors 



The simplest setup of point detector, which guarantees stable reconstruction, 
is the sphere T = <9-Br (0) that encloses the region of interest. Then the set 
of data is 

{p 7 (x, t) | x 6 T, t e [0, T}} (T > sufficiently large) . 

Inserting the spherical mean representation (J7|) into integral equation (}3~Tj) 
and performing integration by parts yield 



— [I(t) * t My(t, t')] 

o 

since (j251) holds and 

'R sp (<j))(x,t' 



{ort>0 



Ant' 



Airt' 



t'=0 



Po(x,t) dr 



0. 



J t'=0 



This implies 

Theorem 3. Let (A1)-(A3) be satisfied and M 7 be defined as in l[21\) . The 

spherical projection R sp (<f)) of satisfies 

t 

(32) P 7 (x,t) = J N^t^Rs^^t')^' x^0,t>0, 

o 

where 



N 7 (t,t') 



Ant 



1 Blt,M 



Case of planar detectors 

In order to define a data setup of planar detectors we need 

Definition 1. Let n G S 1 and E(n,t) denote a plane normal to n with 
normal distance t to the origin. Let be integrable. We define the planar 
projecton of (ft by 

R pl (4>)(n,t) := J 0(x')dA 2 (x') neS\t>0. 
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Let T = S 2 . The simplest setup of planar detectors, which guarantees stable 
reconstruction, is given by 

{Rpiip^-, t))(n, Rq) | n G T, t G [0, T]} (T > sufficiently large) . 

According to [3U [7] the identity 

R pl {p {-,s))(n,Ro) = 2R pl {<j)){n,R - s) n G £\ s > 

holds. This identity and integral equation fl3T|) imply 

t 

R pl (p 7 (-,t))(n,R ) = J 2[I(t) * t M 1 (t,s)}R pl (4>)(n,R -s)ds 

o 

Ro 



2 [/(t) * t M 7 (t, R - t')\ R pl {4>) (n, t') dt' . 
Ro-t 

This leads to 

Theorem 4. Lei (A1)-(A3) be satisfied and M 7 be defined as in l[21\) . The 
planar projecton R pi of (ft satisfies 

Ro 

(33) R vl (p J (;t)){n,R )= J jV 7 (t,t / )i2 p ,(0)(n,Odi' neS\t>0 

Ro-t 

where 

(34) iV 7 (t, t') = 2 [1(f) *t M 7 (t, Ro - 0] • 
Case of line detectors 

In this case the data setup is more complicated. Again we start with a 
definition. 

Definition 2. Let n G S 1 and E(n, 0) be defined as in Definition^ For each 
x£E, / n (x) denotes the line passing through x normal to E. For integrable 
<p we define the line integral operator by 

$ n :=L n (0)(x) := J ^(xOdAV) 

VFe define the circular projection of $ n by 

# circ ($ n )(x,t) := y $ n (x') dA^x') xG£,t>0. 

ast(x)nE 
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Actually we are concerned with three inverse problems. The first is con- 
cerned with the estimation of the circular projections R C irc{&n) for each 
nfS 1 , the second is concerned with the estimation of $ n for each n G S 1 and 
the third is concerned with the estimation of from the set {$ n | n G S 1 }. 
The latter problem corresponds to the inverson of the linear Radon trans- 
form and will not be discussed in this paper. For each n G S 1 let e.g. 
T = <9-Br (0) nE(n), which guarantees a stable reconstruction. Consider the 
sets of data 

M(n) = {L n (p 7 (-,t))(x) |x G T, t G [0,T]} (T > sufficiently large) . 

For each fixed n and data set .M(n) we derive an integral equation for 
Rcird^n)- in P] it was shown that 



(35) 



^ (, S ))(x) = J_|| 



1 d [ R circ ($ n )(x,t') 



Vs 2 - t' 2 



Let M 7 (t, s) := I(t) * t M 7 (t, s). From (IBTj) and (13"5]) and integration by parts 
we get 



o o 



1 

2^ 




t' 



dM 7 (t,s) 



Vs 2 - 1' 2 ds 



ds Rcirc 

($ n )(x,t')dt' 



since M~(t, £) = for £ > and 



1 

2^ 



Lo 



-^circ(^n) (x ; £ ) ^ / 
Vs 2 - t' 2 



s=0 



L n (p (-,r))(x) dr 



0. 



s=0 



We infer 

Theorem 5. Lei (A1)-(A3) be satisfied and M 7 fre defined as in |HJ). .For 
n G S* 1 Zei L n and -R C j rc be defined as in Definition® Then 

t 

(36) L n (p 7 (-,t))(x) = y^ 7 (t,t')^r C ($n)(x,t / )dt / xg£, t>0 



kernel 
(37) iV 7 (M') 



1 

2v7 



Vs 2 - £' 2 



dM 7 (t,s) 
ds 



ds 
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5 Appendix 



In the appendix we prove that the frequency power law (TlOT) for 7 G (1,2] 
fails causality, defined as in Section |2j and that the complex attenuation law 
proposed in (JTT1) permits causality Moreover, we prove that the attenuated 
wave equation has a unique Green function that vanish for t < 0. For this 
purpose we need Theorem 4 in [1], which is stated below (cf. Theorem 7.4.3 
and the following remark in [8]). 

Let M:= {z e C\ Im(z) > 0} denote the upper open complex half plane 
and <S'(R) denotes the space of tempered distributions on R with range in C. 

Theorem 6. A distribution f G <S'(R) is causal, i.e. supp(f) C [0, 00), if 
and only if 

(CI) f : R — > R can be extended to a function F : H — > C that is holomor- 
phic. 

(C2) For all fixed rj > and £ G R, F(£ + ir/) considered as a distribution 
with respect to the variable £ is tempered, and F(£ + i 77) converges (in 
the sense of S' ) when 77 — 0. 

( C3) There exists a polynomial P such that 

\F(z)\ < P(\z\) for Im(z)>e>0. 

If all three conditions are satisfied, then F is the Fourier-Laplace transform 
off- 

Remark 6. The definition of the Fourier transform in this paper has a dif- 
ferent sign as in J^j and thus H is the upper half plane and not the lower half 
plane. 

Theorem 7. For a Q > and 7 G (1,2] let a* be defined as in [Tty) . Then 
F~ 1 {e~ a *^ ' x '}(£) is not a causal function. 

Proof. Let x G M 3 be arbitrary but fixed. 

1) Since e~ a *^^ is a rapidly decreasing function, its inverse Fourier trans- 
form is also a rapidly decreasing function and hence Theorem [6] is applicable. 
The holomorphic power function (— iz) 1 (cf. Remark [2]) defined on C~ is 
the unique holomorphic extension of the function u G R 1— > (— iu) 1 G C and 
hence 

F(z) = exp{-a*(;z) |x|} for z G HI 
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is the unique holomorphic extension of u G R i— > exp{— a*(u) |x|} G C 

2) We show that condition (C3) cannot be satisfied for the sequence (z n ) neN 

defined by 

z n := in G H for n G N. 
From — iz n = n and cos(| 7) < for 7 G (1, 2], it follows 



\F(z ri 



exp < — a 



cos(f 7) 



exp < a 



n ' 



cos(| 7)| 



which cannot be bounded by a polynomial P{n). 

Remark 7.1) We note that the following modification of the power law 

a*(u) = 0:0(7) (-iujy + a x (-iw) (7 G (1,2]) , 
which also appears in the literature, does not satisfies causality, since 

n' 



□ 



F Z r , 



exp < ao 



cos(f 7 )| 



ot\ n 



cannot be bounded by a polynomial for 7 > 1 . 

2) Let 7 G (0,1). It is easy to see that \F(z)\ is bounded by a constant, 
since Re(ce*(z)) is always positive and thus \F(z)\ decreases exponentially 
for \z\ — > 00. 



Theorem 8. For ao, r > and 7 G (1,2], let a* be defined as in 01]) . 
Then J r_1 {e _a *^^ x '}(t) is a causal function. 

Proof. Let x G R 3 be arbitrary but fixed. 

1) Since e _a *^-" x ' is a rapidly decreasing function, its inverse Fourier trans- 
form is also a rapidly decreasing function and thus Theorem [6] is applicable. 
Let M_ e := {z G C|Im(^) > — e}. Since at* (2) is holomorphic on M_ e 
for sufficiently small e, this function is the unique holomorphic extension of 
wGl4 a*(w) e C. Thus 

= exp{-a !)i (2;) |x|} for zeM 

is the unique extension of w G I 4 exp{— a*(u) |x|} G C. This proves 

conditions (CI) and (C2) of Theorem [61 

3) Since the inequality in (C3) is equivalent to 

exp{-Re(a*(z))} < P{\z\) for lm{z) > e > , 
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(C3) is satisfied if 
Since 



Re(a*(H)) C [0, oo) . 



Re(a*) = — (wi z 2 + w 2 z\) 
Co 



with 



/(*) 



1 



+ (-iro^)^ 1 
we have to show that 
(38) w 1 > and 



wi + iw 2 



(z 2 > 0) 



for 



z e M. 



W7 2 > 



Let /i = -ir 2, / 2 (z) := \/l + z^~ l and / 3 (z) := ~, so that / = / 3 o f 2 o /j. 
Moreover, let 



Mi 
M 2 
M 3 
Ni 

We see that 



{^i + iz 2 e C\zi < 0,z 2 > 0}, 
{zi +iz 2 e C|zi > 0,z 2 > 0}, 
{zi +i^ 2 e C|«i > 0,^2 < 0}, 

{i z 2 <E C I z 2 > 0} and N 2 := { Zl e C | z x > 0} . 

/ 1 (M 1 )CM 2 , /i^CJVa, fi(M 2 ) C M 3 , 
/ 2 (M 2 )CM 2 , f 2 (N 2 )CN 2 , / 2 (M 3 )CM 3 . 



From 



1 _ - 1 Z 2 

Zi + i 2 2 + ^1 



+ 2f ^ 



we get 

/ 3 (M 2 )CM 3 , f 3 {N 2 )CN 2 , / 3 (M 3 )CM 2 . 
Therefore we end up with 

/(Mi)CM 3 , f(N 1 )CN 2 , /(M 2 )CM 2 , 

since f = f-s f 2 fi- In other words condition fl38|) is satisfied and hence 
condition (C3) holds. This concludes the proof. □ 

Theorem 9. For aj , t o > and 7 G (1,2], Zet a* 6e defined as in / TJ7]) . 
T/ien equation ITi~b}) with /(x, t) = <5(x) 5(t) and 

(39) P T U<o = 

/ias a unique solution 



p 7 (x,t) = J 7 1 



( exp{ifc(a;) |x|} 
\ 47r|x| 



with k(u) := i a* (a;) + 



Co' 



c 
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Proof. Since p 7 is the Green function, we write G 7 instead of p 1 . The 
Helmholtz equation of f[T5"j) reads as follows 



-2/- , _ <K X ) 



V 2 G^ + fc G- 



\J2n 

and has the only two solutions 



G~(X, = . : : S = ±1 . 

7 v 7 ^ 4tt|x| 

If condition ( 139|) holds, then according to Theorem [6j we have for z = Z\ + 
i 2 2 G EI and x^O: 

(40) e -«R«(«.W)[x| e --«J=f = | e »ifcW|x[| < p(| z |) ZlG M, z 2 >0 
for some polynomial P. From Theorem [S] and Theorem [HI we get 

e -Rfi(a.(*))|x| < Z2 > e > 



for some polynomial Q, i.e. f )40p holds with P := Q for s = 1. If s = — 1, 
then the left hand side of (K0|) grows exponentially, since Re(a*(w)) > is 
increasing, i.e. condition ([39]) cannot be satisfied for s — — 1. Hence there 
exists a unique solution. □ 
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